PhD Chapter 2

Ecological Status calculation and analyses

This series of files compile analyses done during Chapter 2.

All analyses have been done with R 4.0.4.

Click on the table of contents in the left margin to assess a specific analysis.
Click on a figure to zoom it

⏪ | 🏠 | 📄

1. Indicator results

1.1. Scatterplots

Here are the plots presenting the value of each indicator at shallow and deep stations.

1.2. Ecological Quality Status

When relevant, we calculated an Ecological Quality Ratio (EQR) as established by the WFD and MSFD (which varies between 0 and 1). This ratio is calculated with the following equation:

\[ EQR = \frac{V_{ind} - R_{bad}}{R_{good} - R_{bad}} \]

  • \(V_{ind}\) is the value of an indicator at a certain location
  • \(R_{bad}\) is the reference value for a “bad” status
  • \(R_{good}\) is the reference value for a “good” status

This ratio is then classed into Ecological Quality Status (EQS) categories, where reference values and limits for class transitions are specific to each indicator. Five classes are typically described:

  • bad (red #FF0000)
  • poor (orange #FFA500)
  • moderate (yellow #EEEE00)
  • good (green #228B22)
  • high (blue #0000EE)

We calculated this ratio using different indicators, in order to compare their efficiency and relevance.

AMBI

We defined class thresholds using the methods from Borja et al. (2000) and Muxika et al. (2005).

M-AMBI

We defined class thresholds using the method from Muxika et al. (2007).

BENTIX

We defined class thresholds using the method from Simboura & Zenetos (2002).

BOPA

We defined class thresholds using the method from Dauvin & Ruellet (2007).

2. Non-metric Multidimensional Scaling

Species densities have been transformed with a (log+1) operation. Green points are stations within the bay, red points are within the archipelago.

Total density

Total biomass

W-Statistic

Specific richness

Shannon index

Margalef index

Simpson index

Pielou evenness

Taxonomic diversity

Functional richness

Functional evenness

Functional divergence

AMBI

M-AMBI

BENTIX

BOPA

3. Statistical relationships

3.1. Between indicators

In this section, we study the statistical relationships between indicators calculated above.

3.1.1. Covariation

Several types of models were considered to explore relationships: linear, quadratic, exponential and logarithmic. The model with the highest \(R^{2}\) is presented on each plot.

⚠️ Only linear models were implemented for now, as there are some bugs with the calculation of the others.

Total density

Total biomass

W-Statistic

Specific richness

Shannon index

Margalef index

Simpson index

Pielou evenness

Taxonomic diversity

Functional richness

Functional evenness

Functional divergence

AMBI

M-AMBI

BENTIX

BOPA

3.1.2. Correlation

Correlations have been calculated with Spearman’s rank coefficients.

Correlation coefficients between indicators for shallow stations
  N B W S H margalef lambda J delta FR FE FD AMBI M_AMBI BENTIX BOPA
N 1 0.364 -0.359 0.769 0.053 0.252 -0.203 -0.663 -0.44 0.802 -0.327 0.366 0.148 0.373 0.046 0.252
B 0.364 1 0.349 0.426 0.233 0.356 0.15 -0.109 0.082 0.504 0.258 0.294 -0.415 0.483 0.226 0.337
W -0.359 0.349 1 0.131 0.623 0.526 0.675 0.59 0.707 -0.011 0.673 0.112 -0.125 0.368 -0.039 0.147
S 0.769 0.426 0.131 1 0.58 0.759 0.236 -0.333 -0.05 0.873 0.034 0.407 0.111 0.798 0.021 0.445
H 0.053 0.233 0.623 0.58 1 0.811 0.889 0.461 0.59 0.381 0.577 0.215 -0.002 0.778 -0.104 0.298
margalef 0.252 0.356 0.526 0.759 0.811 1 0.613 0.222 0.476 0.583 0.412 0.253 -0.083 0.855 0.062 0.411
lambda -0.203 0.15 0.675 0.236 0.889 0.613 1 0.701 0.754 0.085 0.651 0.078 -0.023 0.502 -0.179 0.083
J -0.663 -0.109 0.59 -0.333 0.461 0.222 0.701 1 0.862 -0.411 0.544 -0.151 -0.296 0.055 0.115 -0.158
delta -0.44 0.082 0.707 -0.05 0.59 0.476 0.754 0.862 1 -0.195 0.514 -0.109 -0.279 0.263 0.065 -0.099
FR 0.802 0.504 -0.011 0.873 0.381 0.583 0.085 -0.411 -0.195 1 -0.066 0.314 0.098 0.636 0.001 0.368
FE -0.327 0.258 0.673 0.034 0.577 0.412 0.651 0.544 0.514 -0.066 1 0.343 -0.223 0.395 0.01 0.207
FD 0.366 0.294 0.112 0.407 0.215 0.253 0.078 -0.151 -0.109 0.314 0.343 1 -0.206 0.431 0.353 0.07
AMBI 0.148 -0.415 -0.125 0.111 -0.002 -0.083 -0.023 -0.296 -0.279 0.098 -0.223 -0.206 1 -0.318 -0.781 -0.184
M_AMBI 0.373 0.483 0.368 0.798 0.778 0.855 0.502 0.055 0.263 0.636 0.395 0.431 -0.318 1 0.253 0.53
BENTIX 0.046 0.226 -0.039 0.021 -0.104 0.062 -0.179 0.115 0.065 0.001 0.01 0.353 -0.781 0.253 1 0.146
BOPA 0.252 0.337 0.147 0.445 0.298 0.411 0.083 -0.158 -0.099 0.368 0.207 0.07 -0.184 0.53 0.146 1
p-values of correlation test between indicators for shallow stations
  N B W S H margalef lambda J delta FR FE FD AMBI M_AMBI BENTIX BOPA
N 0 0.06769 0.07166 4.428e-06 0.7984 0.2149 0.3197 0.0002201 0.02459 8.389e-07 0.1027 0.06626 0.4694 0.06021 0.8253 0.2149
B 0.06769 8.884e-08 0.08109 0.0301 0.2512 0.07501 0.4626 0.597 0.6882 0.008604 0.2029 0.1453 0.03617 0.01241 0.2666 0.09228
W 0.07166 0.08109 8.884e-08 0.5234 0.0008811 0.006435 0.000222 0.001505 8.24e-05 0.959 0.0001671 0.5866 0.5398 0.06469 0.8513 0.4738
S 4.428e-06 0.0301 0.5234 2.315e-189 0.001917 7.044e-06 0.2459 0.09682 0.8101 6.131e-09 0.8675 0.0389 0.5902 1.037e-06 0.9202 0.02268
H 0.7984 0.2512 0.0008811 0.001917 8.884e-08 2.29e-06 1.665e-06 0.01789 0.001851 0.05487 0.002051 0.2906 0.9947 2.858e-06 0.6143 0.1388
margalef 0.2149 0.07501 0.006435 7.044e-06 2.29e-06 8.884e-08 0.001103 0.2759 0.01484 0.00176 0.03665 0.2123 0.6857 2.628e-08 0.7638 0.03696
lambda 0.3197 0.4626 0.000222 0.2459 1.665e-06 0.001103 8.884e-08 6.646e-05 1.509e-05 0.6792 0.000319 0.7063 0.9122 0.008977 0.3818 0.6856
J 0.0002201 0.597 0.001505 0.09682 0.01789 0.2759 6.646e-05 0 1.484e-08 0.03721 0.004045 0.4626 0.1419 0.7906 0.5771 0.44
delta 0.02459 0.6882 8.24e-05 0.8101 0.001851 0.01484 1.509e-05 1.484e-08 8.884e-08 0.34 0.007193 0.5947 0.1677 0.1943 0.7511 0.6316
FR 8.389e-07 0.008604 0.959 6.131e-09 0.05487 0.00176 0.6792 0.03721 0.34 0 0.7474 0.1183 0.6346 0.0004733 0.996 0.06463
FE 0.1027 0.2029 0.0001671 0.8675 0.002051 0.03665 0.000319 0.004045 0.007193 0.7474 0 0.08591 0.2736 0.0456 0.9616 0.3113
FD 0.06626 0.1453 0.5866 0.0389 0.2906 0.2123 0.7063 0.4626 0.5947 0.1183 0.08591 0 0.313 0.02783 0.0765 0.7352
AMBI 0.4694 0.03617 0.5398 0.5902 0.9947 0.6857 0.9122 0.1419 0.1677 0.6346 0.2736 0.313 8.884e-08 0.113 2.494e-06 0.3682
M_AMBI 0.06021 0.01241 0.06469 1.037e-06 2.858e-06 2.628e-08 0.008977 0.7906 0.1943 0.0004733 0.0456 0.02783 0.113 0 0.2133 0.005361
BENTIX 0.8253 0.2666 0.8513 0.9202 0.6143 0.7638 0.3818 0.5771 0.7511 0.996 0.9616 0.0765 2.494e-06 0.2133 0 0.4759
BOPA 0.2149 0.09228 0.4738 0.02268 0.1388 0.03696 0.6856 0.44 0.6316 0.06463 0.3113 0.7352 0.3682 0.005361 0.4759 0

Correlation coefficients between indicators for deep stations
  N B W S H margalef lambda J delta FR FE FD AMBI M_AMBI BENTIX BOPA
N 1 -0.037 -0.306 0.578 -0.044 0.207 -0.226 -0.665 -0.39 0.347 -0.547 0.109 0.156 0.171 -0.102 0.022
B -0.037 1 0.35 0.009 -0.116 -0.002 -0.134 -0.128 -0.133 0.188 0.167 -0.09 -0.119 -0.006 0.077 0.163
W -0.306 0.35 1 0.372 0.747 0.609 0.754 0.627 0.694 0.319 0.419 -0.323 -0.289 0.635 -0.025 -0.074
S 0.578 0.009 0.372 1 0.672 0.902 0.466 -0.057 0.28 0.711 -0.193 -0.267 -0.097 0.792 -0.06 -0.043
H -0.044 -0.116 0.747 0.672 1 0.855 0.956 0.635 0.807 0.462 0.205 -0.386 -0.251 0.873 -0.049 -0.154
margalef 0.207 -0.002 0.609 0.902 0.855 1 0.704 0.288 0.566 0.666 0.064 -0.371 -0.229 0.889 -0.03 -0.078
lambda -0.226 -0.134 0.754 0.466 0.956 0.704 1 0.794 0.893 0.326 0.308 -0.41 -0.284 0.76 0.001 -0.195
J -0.665 -0.128 0.627 -0.057 0.635 0.288 0.794 1 0.878 -0.062 0.586 -0.386 -0.313 0.385 0.075 -0.177
delta -0.39 -0.133 0.694 0.28 0.807 0.566 0.893 0.878 1 0.21 0.432 -0.498 -0.297 0.603 0.019 -0.224
FR 0.347 0.188 0.319 0.711 0.462 0.666 0.326 -0.062 0.21 1 -0.034 -0.281 -0.14 0.575 -0.075 0.027
FE -0.547 0.167 0.419 -0.193 0.205 0.064 0.308 0.586 0.432 -0.034 1 -0.118 -0.199 0.093 -0.097 -0.134
FD 0.109 -0.09 -0.323 -0.267 -0.386 -0.371 -0.41 -0.386 -0.498 -0.281 -0.118 1 0.323 -0.396 -0.237 0.081
AMBI 0.156 -0.119 -0.289 -0.097 -0.251 -0.229 -0.284 -0.313 -0.297 -0.14 -0.199 0.323 1 -0.52 -0.697 0.115
M_AMBI 0.171 -0.006 0.635 0.792 0.873 0.889 0.76 0.385 0.603 0.575 0.093 -0.396 -0.52 1 0.174 -0.097
BENTIX -0.102 0.077 -0.025 -0.06 -0.049 -0.03 0.001 0.075 0.019 -0.075 -0.097 -0.237 -0.697 0.174 1 -0.149
BOPA 0.022 0.163 -0.074 -0.043 -0.154 -0.078 -0.195 -0.177 -0.224 0.027 -0.134 0.081 0.115 -0.097 -0.149 1
p-values of correlation test between indicators for deep stations
  N B W S H margalef lambda J delta FR FE FD AMBI M_AMBI BENTIX BOPA
N 0 0.7422 0.005161 1.332e-08 0.6963 0.06213 0.04084 9.568e-12 0.0002929 0.001393 1.058e-07 0.3297 0.1605 0.1252 0.3631 0.8463
B 0.7422 0 0.001335 0.9369 0.2967 0.9868 0.2305 0.2524 0.2347 0.09033 0.1325 0.4203 0.2875 0.9597 0.4911 0.1432
W 0.005161 0.001335 0 0.0005736 0 1.266e-09 0 0 0 0.003636 0.000103 0.003214 0.008567 1.538e-10 0.8235 0.5089
S 1.332e-08 0.9369 0.0005736 0 4.887e-12 6.707e-31 1.024e-05 0.6104 0.01092 6.969e-14 0.08275 0.01521 0.3857 7.871e-19 0.5898 0.7001
H 0.6963 0.2967 0 4.887e-12 0 1.464e-24 0 0 0 1.565e-05 0.06474 0.0003808 0.02266 1.193e-26 0.661 0.167
margalef 0.06213 0.9868 1.266e-09 6.707e-31 1.464e-24 0 1.684e-13 0.008584 3.017e-08 8.921e-12 0.5687 0.0006022 0.03854 7.53e-29 0.7859 0.4853
lambda 0.04084 0.2305 0 1.024e-05 0 1.684e-13 0 0 0 0.002896 0.005111 0.0001512 0.009821 1.166e-16 0.9929 0.07979
J 9.568e-12 0.2524 0 0.6104 0 0.008584 0 0 0 0.5822 1.402e-08 0.0003808 0.004198 0.0003498 0.5017 0.1119
delta 0.0002929 0.2347 0 0.01092 0 3.017e-08 0 0 0 0.05853 5.912e-05 2.681e-06 0.006821 2.002e-09 0.8653 0.04351
FR 0.001393 0.09033 0.003636 6.969e-14 1.565e-05 8.921e-12 0.002896 0.5822 0.05853 0 0.7584 0.01075 0.209 1.572e-08 0.5053 0.8123
FE 1.058e-07 0.1325 0.000103 0.08275 0.06474 0.5687 0.005111 1.402e-08 5.912e-05 0.7584 0 0.2922 0.07322 0.4079 0.3841 0.2302
FD 0.3297 0.4203 0.003214 0.01521 0.0003808 0.0006022 0.0001512 0.0003808 2.681e-06 0.01075 0.2922 0 0.00305 0.0002311 0.03228 0.469
AMBI 0.1605 0.2875 0.008567 0.3857 0.02266 0.03854 0.009821 0.004198 0.006821 0.209 0.07322 0.00305 0 5.386e-07 3.684e-13 0.3041
M_AMBI 0.1252 0.9597 1.538e-10 7.871e-19 1.193e-26 7.53e-29 1.166e-16 0.0003498 2.002e-09 1.572e-08 0.4079 0.0002311 5.386e-07 0 0.1174 0.3835
BENTIX 0.3631 0.4911 0.8235 0.5898 0.661 0.7859 0.9929 0.5017 0.8653 0.5053 0.3841 0.03228 3.684e-13 0.1174 0 0.1817
BOPA 0.8463 0.1432 0.5089 0.7001 0.167 0.4853 0.07979 0.1119 0.04351 0.8123 0.2302 0.469 0.3041 0.3835 0.1817 0

3.2. Between indicators and abiotic parameters

In this section, we study the statistical relationships between indicators calculated above and different abiotic parameters, in order to understand how well they can be used to detect perturbations.

3.2.1. Covariation

Several types of models were considered to explore relationships: linear, quadratic, exponential and logarithmic. The model with the highest \(R^{2}\) is presented on each plot.

⚠️ Only linear models were implemented for now, as there are some bugs with the calculation of the others.

Total density

Total biomass

W-Statistic

Specific richness

Shannon index

Margalef index

Simpson index

Pielou evenness

Taxonomic diversity

Functional richness

Functional evenness

Functional divergence

AMBI

M-AMBI

BENTIX

BOPA

3.2.2. Correlation

Correlations have been calculated with Spearman’s rank coefficients.

Correlation coefficients between indicators and habitat parameters for shallow stations
  N B W S H margalef lambda J delta FR FE FD AMBI M_AMBI BENTIX BOPA
depth -0.183 0.151 0.361 0.116 0.345 0.252 0.206 0.375 0.361 -0.112 0.21 0.196 -0.324 0.325 0.47 0.022
om -0.073 0.089 0.089 0.107 0.089 0.187 -0.098 0.013 0.022 -0.036 -0.004 0.092 -0.432 0.311 0.449 0.347
gravel 0.096 0.055 -0.026 -0.024 0.071 -0.005 0.196 0.143 0.155 0.223 0.105 -0.194 0.319 -0.151 -0.365 -0.052
sand 0.245 -0.092 -0.016 0.229 0.104 0.082 0.095 -0.118 0.009 0.234 -0.108 0.252 0.471 -0.051 -0.332 -0.273
silt -0.181 0.141 0.041 -0.136 -0.071 -0.03 -0.117 0.049 -0.027 -0.182 0.105 -0.136 -0.474 0.132 0.352 0.252
clay -0.179 0.045 0.021 -0.202 0.048 -0.181 0.168 0.213 0.046 -0.066 0.085 -0.26 -0.119 -0.122 0.008 -0.023
arsenic -0.456 -0.423 -0.1 -0.474 -0.113 -0.182 0.038 0.283 0.225 -0.428 -0.023 -0.595 -0.014 -0.336 -0.124 0.015
cadmium -0.076 -0.416 -0.401 -0.196 -0.182 -0.152 -0.174 -0.066 -0.059 -0.146 -0.301 -0.156 -0.171 -0.11 -0.013 -0.077
chromium -0.522 -0.591 0.022 -0.359 0.033 -0.053 0.06 0.287 0.225 -0.498 0.031 -0.351 0.006 -0.141 -0.019 -0.017
copper -0.546 -0.506 0.138 -0.327 0.05 -0.018 0.068 0.315 0.23 -0.426 0.059 -0.372 0.092 -0.174 0.006 0.011
iron -0.49 -0.391 0.145 -0.309 0.054 -0.006 0.075 0.269 0.262 -0.474 0.083 -0.287 -0.061 -0.103 0.063 -0.037
manganese -0.306 -0.531 -0.109 -0.276 -0.047 -0.073 -0.017 0.092 0.099 -0.459 -0.037 -0.269 0.016 -0.12 0.028 0.041
mercury -0.517 -0.355 0.063 -0.392 0.039 -0.113 0.109 0.423 0.26 -0.499 0.132 -0.379 -0.112 -0.119 0.056 0.021
lead -0.549 -0.499 0.067 -0.381 0.022 -0.034 0.086 0.365 0.237 -0.42 0.062 -0.4 -0.013 -0.174 0.037 0.148
zinc -0.515 -0.613 0.029 -0.343 0.058 -0.037 0.097 0.32 0.259 -0.465 -0.011 -0.411 0.063 -0.172 -0.031 -0.032
p-values of correlation test between indicators and habitat parameters for shallow stations
  N B W S H margalef lambda J delta FR FE FD AMBI M_AMBI BENTIX BOPA
depth 0.3718 0.461 0.06988 0.571 0.08426 0.2148 0.3138 0.05915 0.06988 0.5853 0.3039 0.3373 0.1065 0.1047 0.01535 0.916
om 0.7237 0.6661 0.6636 0.6042 0.6636 0.3587 0.6323 0.9484 0.9175 0.8605 0.9841 0.6537 0.02864 0.1222 0.02133 0.08199
gravel 0.6414 0.791 0.8982 0.9083 0.7293 0.9803 0.3377 0.4863 0.4493 0.2741 0.6109 0.3434 0.1122 0.461 0.06668 0.7995
sand 0.228 0.6549 0.9392 0.2613 0.6145 0.6902 0.6429 0.5648 0.9643 0.249 0.5987 0.2145 0.01519 0.8035 0.09799 0.1778
silt 0.3755 0.4914 0.8422 0.5064 0.7286 0.884 0.5683 0.8112 0.8958 0.3732 0.6092 0.5079 0.01437 0.5192 0.07773 0.2134
clay 0.3828 0.8286 0.918 0.3216 0.8163 0.3756 0.413 0.2965 0.8224 0.7472 0.68 0.1992 0.5637 0.5511 0.9705 0.9109
arsenic 0.01931 0.03116 0.6272 0.01449 0.5839 0.3727 0.8547 0.1608 0.2693 0.02914 0.9105 0.001346 0.9443 0.09292 0.5461 0.9406
cadmium 0.7107 0.03431 0.04234 0.337 0.3728 0.4574 0.3962 0.7489 0.7743 0.4763 0.1345 0.4471 0.4042 0.5918 0.9498 0.7077
chromium 0.006246 0.001481 0.9142 0.07201 0.8722 0.7971 0.7703 0.1547 0.2691 0.009577 0.88 0.07851 0.9762 0.4924 0.926 0.9362
copper 0.003903 0.00909 0.4983 0.1029 0.8098 0.9308 0.7407 0.1171 0.257 0.02999 0.7728 0.06128 0.6539 0.3961 0.9775 0.958
iron 0.01112 0.0489 0.4771 0.1247 0.7917 0.9787 0.7155 0.1831 0.1949 0.01445 0.6878 0.1548 0.7686 0.618 0.7587 0.859
manganese 0.1279 0.005904 0.5945 0.172 0.8202 0.723 0.9334 0.6549 0.6275 0.01828 0.8592 0.1837 0.9388 0.5603 0.8905 0.8406
mercury 0.006883 0.0762 0.7584 0.04779 0.8516 0.5806 0.5945 0.03148 0.1986 0.009534 0.5204 0.0563 0.5852 0.5626 0.7842 0.9176
lead 0.003659 0.009392 0.7461 0.05509 0.9141 0.8709 0.6754 0.06698 0.2429 0.03246 0.7638 0.04293 0.9484 0.3944 0.8577 0.4699
zinc 0.007047 0.0008679 0.8866 0.08579 0.7792 0.8566 0.6358 0.1105 0.2023 0.01656 0.959 0.03695 0.7601 0.4017 0.8813 0.8774

Correlation coefficients between indicators and habitat parameters for deep stations
  N B W S H margalef lambda J delta FR FE FD AMBI M_AMBI BENTIX BOPA
depth -0.17 -0.074 0.3 0.206 0.447 0.309 0.467 0.469 0.534 0.259 0.267 -0.375 -0.166 0.401 0.032 -0.231
om -0.101 -0.037 -0.155 -0.248 -0.13 -0.263 -0.047 0.018 -0.149 -0.206 -0.205 0.009 -0.092 -0.157 0.266 0.203
gravel -0.013 0.205 0.239 0.123 0.076 0.146 0.061 0.058 0.116 0.25 0.095 -0.115 -0.085 0.118 -0.055 -0.023
sand 0.004 -0.015 0.194 0.126 0.179 0.181 0.142 0.111 0.232 0.045 0.178 0.034 0.06 0.15 -0.263 -0.309
silt 0.072 -0.019 -0.219 -0.106 -0.202 -0.187 -0.193 -0.195 -0.286 -0.09 -0.208 0.034 -0.035 -0.139 0.231 0.341
clay -0.019 -0.015 -0.021 -0.004 -0.007 0.011 -0.013 -0.06 -0.044 0.126 -0.092 0.011 -0.034 -0.024 0.041 0.104
arsenic -0.031 -0.079 -0.288 -0.27 -0.285 -0.317 -0.223 -0.169 -0.291 -0.254 -0.204 0.207 -0.048 -0.24 0.225 0.331
cadmium -0.079 -0.077 -0.294 -0.317 -0.289 -0.333 -0.233 -0.154 -0.322 -0.317 -0.184 0.283 0.04 -0.295 0.137 0.276
chromium -0.078 0.034 -0.267 -0.311 -0.329 -0.357 -0.266 -0.174 -0.353 -0.287 -0.221 0.198 -0.026 -0.295 0.235 0.36
copper -0.069 0.006 -0.235 -0.318 -0.285 -0.365 -0.222 -0.169 -0.348 -0.322 -0.252 0.293 -0.06 -0.266 0.254 0.31
iron -0.228 0.072 -0.266 -0.449 -0.363 -0.454 -0.28 -0.094 -0.353 -0.363 -0.189 0.182 0.035 -0.379 0.098 0.328
manganese -0.079 0.049 -0.258 -0.308 -0.313 -0.36 -0.254 -0.163 -0.344 -0.285 -0.275 0.192 -0.03 -0.277 0.196 0.376
mercury 0.019 0.05 -0.203 -0.217 -0.251 -0.278 -0.22 -0.212 -0.329 -0.283 -0.267 0.287 -0.027 -0.205 0.231 0.297
lead -0.041 -0.024 -0.278 -0.301 -0.306 -0.347 -0.256 -0.195 -0.343 -0.273 -0.238 0.283 0.01 -0.282 0.193 0.334
zinc -0.069 -0.053 -0.286 -0.342 -0.308 -0.379 -0.241 -0.182 -0.349 -0.326 -0.218 0.335 0.004 -0.306 0.18 0.303
p-values of correlation test between indicators and habitat parameters for deep stations
  N B W S H margalef lambda J delta FR FE FD AMBI M_AMBI BENTIX BOPA
depth 0.1261 0.511 0.006261 0.06316 2.597e-05 0.004791 9.909e-06 8.77e-06 2.435e-07 0.01902 0.0154 0.0005261 0.1352 0.0001889 0.7763 0.03669
om 0.3675 0.7435 0.1635 0.02459 0.2456 0.01699 0.6754 0.8707 0.1826 0.0639 0.06448 0.9389 0.4133 0.158 0.01573 0.06805
gravel 0.9053 0.06425 0.03055 0.2716 0.4958 0.1918 0.5857 0.6039 0.2993 0.02327 0.3945 0.3057 0.4482 0.2916 0.6212 0.8381
sand 0.9702 0.8968 0.08083 0.2612 0.1078 0.1033 0.2026 0.3188 0.03598 0.69 0.1086 0.7624 0.5932 0.1775 0.01718 0.004724
silt 0.5178 0.8642 0.04799 0.3426 0.06857 0.09176 0.08247 0.07956 0.009136 0.4216 0.06117 0.7603 0.7533 0.2128 0.03641 0.001735
clay 0.868 0.8947 0.8486 0.9691 0.9523 0.9206 0.9067 0.5947 0.6935 0.2609 0.411 0.922 0.7611 0.833 0.714 0.3521
arsenic 0.7815 0.4799 0.008696 0.01399 0.009455 0.003656 0.04401 0.1283 0.007928 0.02149 0.06595 0.06212 0.6656 0.02989 0.04228 0.002378
cadmium 0.4801 0.4899 0.00741 0.003699 0.008482 0.002217 0.03551 0.1676 0.003198 0.003684 0.09801 0.01005 0.7245 0.007097 0.2204 0.01204
chromium 0.4853 0.7593 0.01533 0.004395 0.002539 0.0009817 0.01583 0.1185 0.00114 0.008854 0.04601 0.07398 0.8178 0.007079 0.03379 0.0008835
copper 0.5384 0.9568 0.03393 0.00356 0.009581 0.0007434 0.04513 0.1294 0.001377 0.003164 0.02234 0.007541 0.5946 0.01584 0.02136 0.004539
iron 0.0391 0.5209 0.01611 2.357e-05 0.0008703 1.821e-05 0.01105 0.3981 0.001222 0.0008618 0.08931 0.1023 0.7543 0.0004477 0.3805 0.002633
manganese 0.4802 0.6632 0.01965 0.004816 0.004388 0.0008973 0.0214 0.1436 0.001634 0.009644 0.01279 0.08341 0.7877 0.01168 0.07685 0.0004928
mercury 0.8681 0.6525 0.06723 0.0502 0.02294 0.01157 0.0472 0.05604 0.002681 0.01032 0.01564 0.009027 0.8123 0.06473 0.03687 0.006739
lead 0.7116 0.8275 0.01159 0.005972 0.005188 0.001414 0.0201 0.07908 0.001615 0.01299 0.03097 0.01008 0.9308 0.01013 0.08155 0.002162
zinc 0.5373 0.6367 0.009233 0.001654 0.00481 0.0004453 0.02904 0.1017 0.001332 0.002824 0.04918 0.002117 0.9743 0.005254 0.1064 0.005722

4. Indicator robustness

In this section, we are calculating values of the indicators for shallow and deep stations using a bootstrap method, so that we have an idea of the robustness of each measure.

Bootstrap results for shallow stations
  True mean Bootstrap Mean bias Boostrap 95% CI
S 9.192 -0.03569 [9.1708;9.2852]
N 138.7 -0.3769 [136.8829;141.2555]
B 7.352 -0.05717 [7.1399;7.6789]
W 0.0109 -0.0164 [0.0269;0.0277]
H 1.353 -0.006664 [1.3537;1.3663]
margalef 1.926 -0.01379 [1.9302;1.949]
lambda 0.6202 -0.002792 [0.6205;0.6255]
J 0.6564 -0.00366 [0.6572;0.663]
delta 51.66 -0.3568 [51.7936;52.2472]
FR 23.35 -3.171 [26.111;26.9268]
FE 0.5542 0.002318 [0.5495;0.5543]
FD 0.7656 -0.007396 [0.7701;0.7759]
Bootstrap results for deep stations
  True mean Bootstrap Mean bias Boostrap 95% CI
S 13.99 -0.009295 [13.964;14.0302]
N 89.13 -0.1556 [88.8068;89.7726]
B 8.72 0.0243 [8.5538;8.838]
W 0.02522 -0.007476 [0.0325;0.0329]
H 1.952 0.0002019 [1.9489;1.9551]
margalef 3.046 -7.372e-05 [3.0402;3.052]
lambda 0.7699 0.0001827 [0.7687;0.7707]
J 0.7601 0.0003165 [0.7588;0.7608]
delta 63.48 0.01433 [63.386;63.5518]
FR 31.76 -7.59 [38.827;39.8772]
FE 0.6324 -0.001514 [0.6331;0.6347]
FD 0.8282 0.01091 [0.8165;0.8181]

🔝